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Abstract. The reconnection process in the dynamics of cubic nontwist maps, introduced in [3], is 
studied. The present paper extends the work presented in [8]. As in that work, in order to describe 
the route to reconnection of the involved Poincare-Birkhoff chains or dimerised chains we investigate 
an approximate interpolating Hamiltonian of the map under study revealing again that the scenario 
of reconnection of cubic nontwist maps is different from that occurring in the dynamics of quadratic 
nontwist maps. 
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1. Introduction 

Nontwist maps arise naturally in the study of Hamiltonian systems, because they are models for 
Poincare maps associated to sections in an energy manifold of an iso-energetically degenerate two degree 
of freedom Hamiltonian system |21 1141 IT5] . in transport problems in plasma physics, accelerator physics 
and in other areas. Transport problems in plasma physics can be modelled by an area preserving map 
where the twist condition fails [7|. Applications of nontwist maps in accelerator physics can be found 
in PJ. In 0] are studied the quadratic nontwist standard-like maps both from theoretical and numerical 
point of view. Nontwist standard-like maps exhibit both time-reversal and spatial symmetry being 
observed the appearance of the meanders. Meanders are invariant circles that exhibit foldings in such 
a way that they are not graphs of functions. During the last decade numerical and theoretical studies 
of quadratic non-twist maps [3], EJ, Efl, revealed a global bifurcation (called reconnection) of the 
invariant manifolds of two distinct regular hyperbolic periodic orbits having the same rotation number. 
At the threshold of reconnection the involved hyperbolic orbits are connected by a common arc of their 
invariant manifolds. The physical model of reconnection is met in Tokamaks [13] which are experimental 
machines for achievement of controlled thermonuclear fusion reactions. 

For a rigorous analysis of local and global bifurcations occurring in a family of area preserving maps 
defined on an annulus T x [a, b] (T denotes the circle identified with [0, 2tt)) one derives an approximate 
interpolating Hamiltonian of the map under study 0- 

The present work deals with reconnection in the cubic nontwist area preserving diffeomorphism of the 
annulus TxI,/:(i,i/)h (x', y'): 

, . x' = x + F(a, b; y') (mod 2n) 

y' = y + k sin x 

where the rotation number function F is a cubic map depending on two parameters a > 0, b G R, b ^ 0, 
i.e. F(a, b;y) — y — ay 2 + by 3 and k > is a perturbation parameter. 

We recall that an area preserving diffeomorphism g : T x R — > T x R, g : (x, y) i— > (a/, y') is a twist map 
if d y x' 7^ (d y denotes the partial derivative with respect to y) . Twist property is a basic assumption of 
KAM theorem, as well as for the Aubry-Mather theory pQ. The map (GJ is a non-twist map because it 
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violates the twist condition. Our purpose is to study its dynamics as well as the route to reconnection 
in the case when the shape parameter a and the perturbation parameter k are fixed and the other shape 
parameter 6 varies on the real line or on an interval. 

2. Properties of cubic nontwist map 

First we recall some properties of the map under study jS] . The motion in the unperturbed map , 
i.e the map corresponding to k = 0: 

x' = x + y — ay 2 + by 3 (mod 2ir) 
V = V 

occurs along the circle y — est. The rotation number of an orbit starting at (x, y) is: 

X n — X 



(3) 



lim 

n — *oo 



2nir 



F(a,b;y)/(2n), 



where (X n , Y n ) is the orbit of the point (X, Y) = (x, y) under the lift of the map (the map defined on M 2 
having the same expression, without modulo 2ir for the first component). The map (J2J violates the twist 
condition for the parameter values (a, b) such that a 2 — 36 > 0, along the circles: 



(4) 



Ci: y = 



a + \fa? 



36 



C* 2 



y = 



a — \J o? — 36 



36 ' * " 36 
These circles are called twistless or shearrless circles. At the same time along the circle C\ the rotation 
number has a global minimum, 

pi _ -2a 3 -3a 2 ^a' 2 -3b+9ab+(a 2 +6b)^a' ! -3b 
r ICi — 5 4Trb 2 ^ , 

while along C*2, a global maximum, 



F 



c> 



- 2a 3 +3a 2 Va 2 -3b+9afc- (a 2 +6b) \/a 2 -3b 



5iiTb 2 

Let us denote by u) m ,u>M' 



(5) 



a + ya 



3b 



Va 2 - 36 
~^6 ' 



the points of minimum, respectively maximum, for the rotation number function F. For y G (— oo,cjm) 
the unperturbed map has positive twist (the rotation number function is increasing), for y £ (cuMi^m) 
has a negative twist (the rotation number function is decreasing), while for y S (w m ,+oo) has again a 
positive twist (Fig^l. 
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Figure 1. The graph of the rotation number function F(a, 6; y) for a = 2.5 and 6 = 1.26 
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The orbits lying on the circles y — yo with F(a,b;yo)/2n a rational number in lowest terms, p/q, 
are periodic orbits. If such a periodic orbit lies in a region of monotone twist property of the map, 
after a slight perturbation, it gives rise generically to at least two periodic orbits of the same rotation 
number, one elliptic and the second regular hyperbolic. Elliptic points are surrounded by invariant circles, 
and hyperbolic points are connected by heteroclinic connections. Such a pair of periodic orbits and the 
associated invariant sets form a Poincare-Birkhoff chain. For y e (com, w m ) can exist three circles y = est 
of the unperturbed map, on which lie periodic orbits of the same rotation number p/q. Our aim is 
to study the bifurcations of the periodic orbits or the invariant manifolds belonging to three distinct 
Poincare-Birkhoff chains created after a slight perturbation, as the shape parameter 6 € R, a 2 — 36 > 
defining the rotation number function F varies. 

3. RECONNECTION SCENARIO 

In order to analyze the changes in the topology of invariant manifolds of the involved p/q- type hyper- 
bolic periodic orbits consider the interpolating Hamiltonian associated to the map F: 

(6) H aAk (x,y) = -y 2 /2 + ay 3 /3- by 4 /i-kcosx 
It defines the vector field 

(7) X H = f — , — \= (y- ay 2 + by 3 , ksinx) 

which is reversible with respect to the involution R(x, y) ~ {—x, y), i.e. R o Xh = —Xh o R. The fixed 
point set, Fix(i?), consists in the lines x = and x = it, called symmetry lines. The equilibrium points of 
Xh lying on the symmetry lines are called symmetric. The Hamiltonian system associated to the vector 
field Xh can display at most three chains: Poincare-Birkhoff chains or dimerised chains. A dimerised 
chain is a structure formed by elliptic points surrounded by homoclinic circles to the corresponding 
hyperbolic points. 

In order to describe the scenario of reconnection and the local bifurcations of the equilibrium points, we 
analyze the position on the symmetry lines of the equilibrium points, their stability type and bifurcations 
occurring as b varies and a, k are fixed in the parametric space (a, b, k), with a, fc > and 6 £ 1. 
Therefore, if b < the Hamiltonian system has six equilibrium points (e stands for elliptic and h for 
hyperbolic): 

Pih(0, 0); P 2e (7T, 0); P 3e (0, a ~ V ^ rb ); 

PAh{-K, a ' V ^ rb ); P 5h (0, °+vCT ) ; p 6e (^ o+vgEH) 

If ^ < b < ^ (Fig[5J the vector field Xh has only two equilibrium points: Pi(0, 0), Piij^-, 0), while for 
a 2 — 46 it has four equilibrium points: Pi(0, 0), P2(tt, 0), ^4(0, B(tt, -^=). In the latter case, a 2 — 46 = 0, 
the two eigenvalues are zero and a bifurcation of equilibrium points occurs. 

In the following we want to describe the local changes in the topology of the invariant manifolds of 
the Hamiltonian system, when a, k > are fixed and 6^0 varies on the real line. We remark that to 
get connected any two neighboring chains when 6 varies, we need to consider both positive and negative 
values of the parameter 6. When the parameter a varies [S], it is sufficient to consider the case a > 0. The 
systems whose phase portraits are illustrated in different figures correspond to a = 1.5 and fc = 0.018. 
Denote by I, II and III the three chains containing the equilibrium points, more precisely, the chain I 
contains the points P\h, Pi e , the chain II contains the points P% e , P^h and III the points P^h, P&e- For 6 
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Figure 2. The existence domain of the equilibrium points 

small enough (b < —2 for example), the all six equilibrium points are born, Fig|3^). The chains I and II are 
two dimerised chains. Between these chains the trajectories of the Hamiltonian vector field are not graphs 
of real functions of x, but they are meanders. Each two neighboring points on the same symmetry line 
have opposite stability type. Increasing further the parameter b, the equilibrium points lying on the same 
symmetry line (such points lie within two different chains) go away. At a critical value, called threshold 
of reconnection, the hyperbolic points of the two neighboring chains get connected by common branches 
of their invariant manifolds, FigEt>). The common branches are W u {P\h) and W s (Pih) respectively 
W s (Pih) and W u (Pih)- In order to get this threshold of reconnection we impose that the hyperbolic 
equilibrium points P\h and P±h to belong to the same energy level set, that is, H a ^,k(.Pih) = Ha.b.k{P^h)- 
This implies that the reconnection surface of the dimerised I and II chains is: 

(8) 6b 2 + a 4 - 6a 2 b + 48b 3 k + \ab\j a? - 46 - a 3 \j a? - 4b = 

Numerically, it leads to the first threshold of reconnection b := b\ rec = —1.9538. 

Increasing b slightly from &i rec , the two dimerised chains become two Poincare-Birkhoff chains, FigO:), 
so the system displays three Poincare-Birkhoff distinct chains, Fig^Ji) and the Poincare-Birkhoff chains 
II and III approach each other. At the threshold of reconnection the hyperbolic points of these two chains 
(II and III) get connected by common branches of invariant manifolds Fig^Jj). The common branches 
are W u {P ih ) and W s (P 5h ) respectively W s {P ih ) and W u (P 5h ). As above, imposing H a ^ b , k (P ih ) = 
H a ,b,k{P5h) we find the surface of reconnection of the chains II and III: 

(9) fc=l^(a 2 - 46)7^46 

For the numerical values, the threshold of reconnection is b := £>2rec = 0.53168, Fig^Js). Continuing to 
increase b beyond &2rec, the Poincare-Birkhoff chains II and III are transformed into two dimerised chains, 
FigEt). Between these chains, the same as above for the dimerised I and II chains, the trajectories of 

2 

the Hamiltonian vector field are meanders. At the value b = the equilibrium points P^,^P^,P^ and 
Pq are reduced to two points ^4(0, -^=), B(-k, -^=), Fig^Ji), which completely disappear for ^- < b < 
(remember that we work only on the nontwist domain, b < Fig^t). 




Remark that we can simultaneously get connected the all three chains I, II and III. Call this triple 
reconnection. Let us describe the triple reconnection scenario. Consider in this case b > 0. For b slightly 
beyond 0, the chains I and II are two dimerised chains while III is a Poincare-Birkhoff chain, FigEJi). 
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Figure 4. Reconnection scenario of the chains II and III. The values of the parameter 
b are: a) b = 0.5 b) b = 0.53168 c) b = 0.54 d) 6 = 0.5625 e) b = 0.6 



Increasing further the parameter b, the equilibrium points of the two dimerised chains lying on the same 
symmetry line go away while the equilibrium points of the Poincare-Birkhoff chain approaches the points 
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of the dimerised chain II. At a critical value, called triple threshold of reconnection, the hyperbolic points 
of the three chains get connected by common branches of their invariant manifolds FigEb). Imposing 
H ai b,k(Pih) = H a: b,k(P4h) = H a ,b,k{P5h)i we g et from © and © the following reconnection curve: 

(10) a 4 - 6a 2 b + 66 2 + a{a 2 - 46) - 46 = 

Solving numerically (|10[1 for a = 1.5 one get 6 := 63^ = 0.5 and from JHJ k := k^ rec = 0.0625. 
Consequently, to get connected the three chains I, II and III, in the case when 6 varies, we have to keep 
the perturbation parameter k at the constant value k^ rec . For 6 > b^ rec the scenario is similar to the case 
b > b 2 rec, FigEt)-e) and Figit)-e) 

4. Conclusions 

In this paper we have extended the studies reported in [Hj on reconnection scenario of a three-parameter 
cubic nontwist map depending on the parameters o, 6 and k. Using an approximate interpolating Hamil- 
tonian of the map we have described the reconnection process of any two neighboring chains in the case 
when the parameters a, k are fixed and 6 varies. By numerically computations we found the exact values 
of the thresholds of reconnection. At the end we presented the triple reconnection of the all three involved 
chains. 
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